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ARTUR AVILA AND SVETLANA JITOMIRSKAYA 

Abstract. We develop a quantitative version of Aubry duality and use it to obtain several sharp 
estimates for the dynamics of Schrodinger cocycles associated to a non-perturbatively small ana- 
lytic potential and Diophantine frequency. In particular, we establish the full version of Eliasson's 
reducibility theory in this regime (our approach actually leads to improvements even in the per- 
turbative regime: we are able to show, for all energies, "almost reducibility" in some band of 
analyticity) . We also prove 1/2-Holder continuity of the integrated density of states. For the 
almost Mathieu operator, our results hold through the entire regime of sub-critical coupling and 
imply also the dry version of the Ten Martini Problem for the concerned parameters. 



1. Introduction 

This work is concerned with quasiperiodic Schrodinger operators H = H\ v a g defined on I (Z) 
(1.1) (Hu) n = u n+ i + u„_i + Xv(6 + na)u n 

where v : R/Z — > R is the potential, A € R is the coupling constant, a € R \ Q is the frequency and 
9 £ R is the phase. The central (and in some sense most important) example is given by the almost 
Mathieu operator, when v(x) — 2cos(27rx). 

Except where otherwise noted, below we assume the frequency a to be Diophantine in the usual 
sense, and v analytic. 

An important feature of quasiperiodic operators is that the family {-ffAu.a.ejAeR undergoes a so 
called metal- insulator transition when |A| is changed from small to large. Loosely speaking, large A 
correspond to the insulator phase, with positive Lyapunov exponents and related localization type 
effects, while small A leads to a metallic phase, with zero Lyapunov exponents and good transport 
properties. Here we are interested in the metallic regime, and therefore in small couplings. 

One should distinguish between two possible regimes of small |A| (similar considerations can be 
applied to the analysis of large coupling). One is perturbative, meaning that the smallness condition 
on |A| depends not only on the potential v, but also on the frequency a: the key resulting limitation 
is that the analysis at a given coupling, however small, has to exclude a positive Lebesgue measure 
set of a. Such exclusions are inherent to the K AM- type methods that have been traditionally used 
in this context. The other, stronger regime, is called non-perturbative, meaning that the smallness 
condition on |A| only depends on the potential, leading to results that hold for almost every a. Let us 
stress that, in certain related contexts (multifrequency modifications of (jl.ll) ). perturbative results 
do fail to extend to the non-perturbative regime, see Remark II. II 

In [E] , Eliasson obtained a very precise description of operators (|1.1[) in the case of small analytic 
potentials in the perturbative regime. He proved fine estimates on all solutions of the eigenvalue 
equation Hu = Eu for E in the spectrum, concluding that most, but not all, are analytic Bloch 
waves, i.e. quasiperiodic and analytic in the hull. In his context, the problem of existence of analytic 
Bloch waves can be restated (and is indeed treated) as a dynamical systems problem, reducibility 
of the associated cocycle (see §1.2| . His method is based on a sophisticated KAM scheme, which 
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avoids the limitations of earlier KAM methods (that go back to the work of Dinaburg-Sinai [DiS 
and that excluded parts of the spectrum from consideration). 

More recently, a less precise analysis of small analytic potentials has been carried out through the 
non-perturbative regime in [J] and [BJlj . One feature of those results is that most of the analysis is 
concerned with the dual model H = H\ v , a .e defined on l 2 (Z) 

(1.2) (Hu) n = XikUn-k + 2 cos(27T0 + na)u n , 

where Vk are the Fourier coefficients of v(x) = ^ i!ke 2n ' lkx . More precisely, localization (pure point 
spectrum with exponentially decaying eigenfunctions) results for the family {i?AD,a,e}eeR are used 
to obtain information on the family {i?Au,a,e}eeffi- At the root of this approach is the classical Aubry 
duality: a Fourier-type transform corresponds localized eigenfunctions Hu = Eu to analytic Bloch 
waves for the equation Hu — Eu. Duality, originally discovered in [AAj has been given several 
rigorous interpretations since the early 80s. Dynamical version of Aubry duality is that 
localization for the dual model 

leads to reducibility for almost every energy |P2j . A more subtle duality statement is that pure 
point spectrum for the dual model allows to conclude purely absolutely continuous spectrum for a.e. 
9 [GJLSj . However, all the duality links that have been established so far were entirely algebraic, 
with no quantitative estimates involved. Thus they could not lead to an understanding of the whole 
spectrum, since they do not acknowledge the set of energies where localization (and reducibility) 
fails. Inability to study this (zero measure yet topologically generic) set of bad energies was a major 
stumbling block in answering various spectral questions, including some long open conjectures. For 
example, a question whether there is a possibility of collapse of potential gaps, whether there is ever 
any singular spectrum, or what is the exact modulus of continuity of the integrated density of states 
all require estimates holding for all energies. 

In this paper we address this issue by developing the first quantitative version of duality, that 
makes it possible to obtain fine dynamical estimates from local information on the dual model. We 
then show that required local properties do hold non-perturbatively for all energies in the spectrum, 
thus allowing us to obtain 

all of Eliasson's results (and more) in the non-perturbative setting. As discussed above, results 
"for almost every energy" can be obtained by considering only the localized solutions for the dual 
model. Our achievement here is in developing the technique that can handle all energies, including 
those for which localization does not hold, and there is no reducibility. For those energies, our 
approach still gives tight dynamical estimates that imply "almost reducibility" : though coordinate 
changes can not trivialize the dynamics (this would be reducibility), they come arbitrarily close to 
it (Theorem 1 1.4[) . Here the notion of closeness is rather strong, involving control in a fixed band of 
analyticity (such strong control is new even in the perturbative regime). 

Besides almost reducibility, the dynamical estimates yield almost immediately a number of corol- 
laries, including sharp estimates such as l/2-H61der continuity of the integrated density of states 
(Theorems 11.21 and II. 6|) . and the non-collapse of spectral gaps for the almost Mathieu operator 
(Theorem 1 1.1[) . Certain other consequences are less immediate and will be reported separately. 

Our almost reducibility result obviously implies that the non-perturbative setting can be reduced 
(via coordinate change) to the perturbative regime^ so the known results in the perturbative theory 
become automatically non-perturbative. This is applied to conclude absolute continuity of spectral 
measures for all phases (Theorems 11.31 and II. 5|) using a result of [E] . While this reduction also gives 
other corollaries, we stress that Theorems 11.11 11.21 and 11.61 are obtained here without using this 
route. 



1 Theorem 1 1.4 1 is actually much more than is needed for mere reduction to the perturbative regime, which in itself 
is a less delicate result (it is presented in the appendix as Theorem I A. 2) . 
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Remark 1.1. In [Ej, Eliasson is actually able to deal with multifrequency potentials. It was shown 
by Bourgain |B2j that Eliasson's results, in the multifrequency case, do not hold non-perturbatively. 

Remark 1.2. A further common advantage of non-perturbative approaches, which we will not pursue 
here, is to extend through much weaker Diophantine conditions than what can be covered by KAM 
based methods, which usually stop working at the Brjuno condition. The possibility of going beyond 
Brjuno is very interesting for certain problems, see for instance [AJlj . 

1.1. The almost Mathieu operator. Before stating precisely the non-perturbative version of 
Eliasson's result, we discuss the almost Mathieu operator. Results are particularly neat for this case 
because the dual model of the almost Mathieu operator is a rescaled almost Mathieu operator (with 
the same frequency, but inverse coupling), thus the non-perturbative analysis extends through the 
whole subcritical regime |A| < 1 (and can not be extended further, as the almost Mathieu operator 
undergoes a phase transition at |A| = 1 | AAj ) . 

In [AJlj . it was proved that the spectrum of the almost Mathieu operator is a Cantor set for 
any q e R \ Q, A ^ 0. This is the Ten Martini Problem of Kac-Simon. There is a much more 
difficult problem, which remains open, known as the dry version of the Ten Martini Problem: it 
asks to show that the spectrum is not only a Cantor set, but that all gaps predicted by the Gap- 
Labelling Theorem |BLT1 [JM] are open (see Section [4~2l for a precise formulation). In contrast with 
the Ten Martini problem that only requires certain property to hold densely in the spectrum, this 
formulation requires handling of all energies with rational rotation numbers, thus doesn't leave any 
room for energy exclusion. The dry Ten Martini Problem enjoyed a significant attention since its 
first formulation in the early 80s. An affirmative answer was obtained for Liouville a [CEYj . as well 
as for a set of (A, a) of positive Lebesgue measure [PI] . Here we are able to deal with almost every 
(A, a). 

Let DC(k, r) be the set of a G K such that a—- > nq~ T , p, q G Z, q ^ 0. The set of Diophantine 
numbers is DC = U k> o. t >oDC(k, r). 

Theorem 1.1. The dry version of the Ten Martini Problem holds for a G DC, A ^ —1,0, 1. 

The work of Puig mentioned above obtains the same result in the perturbative regime a G DC, 
In |A| large (depending on a) and is based on the perturbative results of Eliasson. 

We note that our proof is not dependent on [AJlj , where the issue was in handling the arithmeti- 
cally critical non-Diophantine regime, while here we entirely focus on the Diophantine case. 

As discussed in [AJlj . in order to bridge the gap between generic and full measure one must 
analyze a particularly difficult region of parameters which is still badly understood (and which this 
work does not touch). The case A| = 1 is also very open even in the case when a is the golden mean 
(the problem is that one can not use localization methods). 

We next move our focus to the regularity of the integrated density of states, a recurring theme 
in the analysis of quasiperiodic operators. 

Theorem 1.2. Let a G DC, A =/= —1,0,1. Then the integrated density of states is 1/2-Hdlder 
continuous. 

This estimate is optimal in several ways. First, there are square-root singularities at the bound- 
aries of gaps (e.g., |P2j ). so the modulus of continuity cannot be improved. Also, it is known that 
for a certain non-empty set of a which satisfies good Diophantine properties (but has zero Lebesgue 
measure) and A = 1, the integrated density of states is not Holder f |B3j . Remark after Corollary 
8.6). Finally, for any A ^ and generic a, the integrated density of states is not Holder (this is 
because the Lyapunov exponent is discontinuous at rational a, which easily implies that it is not 
Holder for generic a). 
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Goldstein-Schlag |GSj had previously shown (1/2 — e)-H61dcr for any e, all A ^ — 1, 0, 1 and a full 
Lebesgue measure subset of Diophantine frequencies (they are also able to consider other analytic 
functions, in the regime of positive Lyapunov exponent, at the cost of a worse Holder constant). 
Previously Bourgain |B1| had obtained almost 1/2- Holder continuity in the perturbative regime, 
for Diophantine a and In |A| large (depending on a). More recently Sana Ben Hadj Amor obtained 
l/2-H61der continuity in the perturbative regime of Eliasson |Amj . 

Another consequence of our results is the following. A well known conjecture, dating back to the 
work of Aubry- Andre [AAj , and more recently included in the list of problems of Simon (Problem 6 
of |S2j). asks to show that for < |A| < 1 the spectral measures of the almost Mathieu operator are 
absolutely continuous. This is also tied to another, more general and far-reaching, conjecture, that 
is sometimes attributed to Simon, saying that for almost periodic operators, singular spectrum must 
be phase independent. While a.c. phasc-indcpcndcncc of the spectrum and its a.c.,s.c. and pure 
point components is an almost immediate corollary of ergodicity, almost periodic operators exhibit 
certain phase rigidity, in that the spectrum |ASj and absolutely continuous spectrum [LSj are the 
same for all phases. Although, this is not true for singular continuous and pure point components 
taken individually (there may be a dependence on the arithmetics of the phase |JSj ). the conjecture 
is that combined together they will also demonstrate the phase stability. This question is only 
non-trivial (but highly so) in the regime of zero Lyapunov exponents. We should note, that by the 
time of the inclusion of the above question as Problem 6 in |S2j . it was already solved in case of 
Diophantine frequencies for a.e. phase [J], thus, in the Diophantine regime, the issue was precisely 
the phase rigidity. This is what we address here. 

Theorem 1.3. The spectral measures of the almost Mathieu operator are absolutely continuous for 
a £ DC, < |A| < 1, 9 £ M. 

In the perturbative regime, a £ DC, < |A| < Ao(a), the result for all 9 £ M. follows from the 
work of Eliasson. As was already mentioned, Theorem 1 1.31 is obtained here from Eliasson's result by 
"non-perturbative reduction to the perturbative regime" . 

Remark 1.3. Though not relevant to this work (since the estimates involved in our reduction to 
the perturbative regime are already bound to the usual Diophantine condition), we should point 
out that this approach does limit possible extensions beyond the Brjuno condition (see Remark 
II. 2j) . Such issues are particularly relevant in view of recent progress towards absolutely continuous 
spectrum "from the Liouville side" [ADj . which seems to break down strictly beyond the Brjuno 
condition^, while absolutely continuous spectrum is expected to hold, as described above, without 
any exceptions 

1.2. Almost reducibility. As previously discussed, Eliasson's analysis is based on dynamical sys- 
tems considerations. A cocycle is defined by a pair (a, A) where a£l and A : R/Z — » SL(2,R) is 
analytic. It is viewed as a linear skew-product (x, w) i— » (x + a, A(x) ■ w), x £ R/Z, w £R 2 . We say 
that two analytic cocycles (a,A^), i = 1,2, are analytically conjugate if there exists an analytic 
map B : M/Z -» PSL(2,M) such that 

(1.3) A {2 \x) = B(x + a)A^ix)B(x)- 1 . 

The dynamical properties of cocycles are preserved by conjugacies. We say that a cocycle is reducible 
if it is C"-conjugate to a cocycle of the form (a, A*) where A* is a constant matrix. Eliasson's 
reducibility theory describes the dynamics of (a, A) when a is Diophantine and A is close to a 

2 Namely, the method of |AD| can only cover irrational numbers that are exponentially well approximated by 
rationals. 

After this work was completed, the first author has announced a complete solution to Problem 6 of S2 , which 
does involve the development of a fully non-perturbative approach to absolutely continuous spectrum. 
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constant. He shows that such cocycles are typically (in a measure-theoretic sense) reducible, and 
gives good estimates for the non-reducible ones. The precise closeness quantifier defines the Eliasson's 
perturbative regime. See the appendix for a summary of the results of the theory that are relevant 
to this work. 

Eliasson's perturbative regime is not invariant under conjugacies. A more intrinsic notion, almost 
reducibility, introduced by Avila-Krikorian AK2j^ in the smooth category, captures the properties 
of cocycles "that behave as cocycles in Eliasson's perturbative regime" . The results of this work 
justify introducing the corresponding definition in the analytic category as well. 

Definition 1.1. An analytic cocycle (a, A) is C"-almost reducible if the closure of its analytic 
conjugacy class contains a constant. 

If (a, A) is C^-almost reducible and a is Diophantine, then (a, A) is analytically conjugate to a 
cocycle in Eliasson's perturbative regime. 

The connection to Schrodinger operators is clear when 



properties of the family of Schrodinger operators {H\ Vt Q .e}eeR are closely related to the dynamics 
of the family of cocycles {(a, S\ Vi e)}e£Wl- 

Theorem 1.4. For a £ DC, v : K./Z — > K analytic, and < |A| < \q{v), the cocycles associated 
with {Hx v ,a,d}eem are almost reducible. If v{x) = 2cos27ra; then Xq = 1. 

Direct application of Eliasson's reducibility theory immediately yields generalized version of The- 
orem |TT3j 

Theorem 1.5. Let a,v,X be as in Theorem \l-4\ Then the singular spectrum is the same for all 
phases 9 £ M (and empty). 

We note that in the non-perturbative regime, the a.e. absolutely continuous spectrum was estab- 
lished in [BJlj . What we address here is the stability of singular spectrum. 

1.3. Further non-perturbative analysis. 

Theorem 1.6. Let a, v, X be as in Theorem \1.4\ Then the integrated density of states is 1/2-Holder. 

The history of this question is discussed after Theorem 1 1.2[ and as for the almost Mathieu case, 
this result is optimal in several ways. While we give a direct simple proof, this theorem can also be 
derived via Theorem 11.41 from a recent perturbative result of Sana Ben Hadj Amor [Amj . 

Our analysis also allows to study a more delicate question: Holder continuity of the individual 
spectral measures (of which the integrated density of states is an average). We can show that for 
all 9 and all localized initial vectors their spectral measures are uniformly l/2-H61der. This result is 
once again optimal, and new even in the perturbative regime. This corollary is a little more involved 
and also needs the introduction of some additional theory. It will be reported separately |AJ2j . 

The dry version of the Ten Martini Problem is specific to the almost Mathieu operator, and does 
not hold for general analytic potentials [DJ] , What can be concluded 

is that for the great majority of v (that is, excluding a set of infinite codimcnsion), and except 
for countably many < |A| < Xq(v), all gaps are open. 



See also the first (preprint) version of | AK1| . 



(1.4) 





In other words, the spectral 
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I. 4. Almost localization and quantitative Aubry duality. Our key problem is thus to show 
almost reducibility of certain cocycles. 

As mentioned before, there is a classic Aubry duality link between localization and reducibility 
(Theorem I2.5[) , whose application is however limited since localization (pure point spectum with 
exponentially decaying eigenfunctions) in general does not hold for every 9 |JS| . This of course fits 
with the fact that reducibility in general does not hold for all cocycles. (Localization for almost 
every 9, proved in [J] and [BJlj . turns out to be enough for many results, but not for the fine ones 
we are interested in here.) 

It is thus natural to devise a weakened notion of localization that could be expected to hold 
for every phase, and to develop ways to link it to almost reducibility. Here we show that this 
approach indeed works. Namely we establish a quantitative version of Aubry duality that links 
local exponential decay of solutions to the dual eigenvalue problem to fine dynamical estimates, 
thus showing that almost reducibility of cocycles associated {i?At),Q,e}eeffi can De concluded from 
a property of the dual model {H\ v , a .e}, which we call "almost localization" (see Definition 13. ip . 
Informally, almost localization gives a precise description of the decay of generalized eigenfunctions 
away from a sparse sequence of resonances, somewhat similar to what is considered in [JLj . but 
more precise than what is obtained in Q] and [BJlj . Refining those results, we indeed show in 
Theorem 13. 21 that almost localization holds for the dual model in the regime relevant for Theorem 

II. 41 Then we proceed to show that for a £ DC, almost localization for the dual model implies almost 
reducibility. This link is quantitative and is significantly more subtle than that between reducibility 
and localization. The key dynamical estimates needed to establish this link, Theorems 13.41 and 13. 5( 
coupled with Theorem 13.21 immediately imply Theorem 1 1 -41 T and all its corollaries). They also imply 
directly Theorems 11.11 11.21 11.61 Those direct implications are not difficult but involve conceptually 
new arguments, presented in Sj4j Particularly the direct proof of Theorem 11.11 is based on linking 
resonant rotation numbers to resonant phases (see Theorem 14. 2\ or Theorem 14.11 for more detail). 
We note that the proof of Theorem 13. 21 (Sj5| is building on well developed localization methods, and 
borrows a number of ingredients from |BJ1] and some shortcuts from [AJ1] . Theorems 13.41 and 13.51 
and their proofs, presented in sections and [HJ are the main novel technical and conceptual 
contributions of this paper, and are not based on previous work. 

Remark 1.4. Previous progress in extending Eliasson's results by coupling non-perturbative and 
perturbative methods was obtained by Avila-Krikorian. In |AKlj , almost every energy results were 
obtained following a rather different technique (reduction to the perturbative regime was achieved 
by renormalization) . More precise (unpublished) results were later obtained (using the connection 
between localization and reducibility in both directions), including a non-perturbative proof of re- 
ducibility under a topological condition (which covered "most" of the spectrum but excluded the 
gap boundaries). 

Let us conclude with a few comments on the analysis of general operators , without restriction 
on the coupling constant!! There has been much recent progress on the description of the part of the 
spectrum E + corresponding to energies with a positive Lyapunov exponent (see [B3| for a lengthy 
account, and |GS| for more recent results). The concept of almost reducibility allows us to determine 
another region of the spectrum E ar which can be thoroughly analyzed (either directly by the methods 
developed here, or by reduction to the perturbative regime). 

There is a number of parallels between our work and previous developments regarding the positive 
Lyapunov exponent regime. 



Earlier, less complete (due to reliance on perturbative techniques), considerations (in the case of smooth poten- 
tials) first appeared in the preprint verison of |AK1) . 
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For instance, the stability of positivity of the Lyapunov exponents (in other words, the property 
that the set of energies corresponding to zero Lyapunov exponent is closed) |BJ2j is paralleled by the 
fact that it is also possible to conclude stability of almost reducibility from Theorem 11.41 utilizing 
an argument of |AK2] . It is precisely the stability of almost reduciblity that makes the concept so 
interesting from a dynamical systems point of view. 

It is interesting to note that almost reducibility may be after all determined by a "mirror" condi- 
tion to positivity of the Lyapunov exponent. Indeed, for energies in the spectrum, almost reducibility 
clearly implies strong vanishing of the Lyaupunov exponent, in the sense that the cocycle may grow 
at most subexponentially in some band. We conjecture that the converse also holds, that is, almost 
reducibility should follow from strong vanishing of the Lyapunov exponent. 

2. Preliminaries 

For a bounded analytic function / defined on a strip {\5sz\ < e} we let ||/|| e = sup| Qz | <£ |/(#)|. If 
/ is a bounded continuous function on R, we let ||/||o = sup^gjj |/(^)|- 

2.1. Cocycles. Let a£l\Q,i£ C°(R/Z, SL(2, C)). We call [a, A) a (complex) cocycle. The 
Lyapunov exponent is given by the formula 

(2.1) L(a,A) = lim - / In ||A n (s)||dx, 

n— >oo n J 

where A n , n £ Z, is defined by (a, A) n — (na, A n ), so that for n > 0, 

(2.2) A n (x) =A(x+ (n- l)a)---A(x). 
It turns out (since irrational rotations are uniquely ergodic), that 

(2.3) L(a,A) = lim sup - In ||A„(a;)||. 

Remark 2.1. By subadditivity, for any compact set K C (R \ Q) X C°(R/Z x SL(2,R)), for every 
5 > there exists Ck,s > such that for every k > 0, 

(2.4) sup sup ln\\A k {x)\\ < C K ,s + k{ sup L(a,A)+5). 

(a,A)€K xGR/Z (a,A)eK 

We say that (a, A) is uniformly hyperbolic if there exists a continuous splitting C 2 = E s (x)®E v (x), 
x e R/Z such that for some C > 0, c> 0, and for every n > 0, • w\\ < C7e _c "||w||, w £ E s (x) 

and || A_„(x) • w\\ < C7e^ cn ||w||, w £ E u (x). In this case, of course L(a,A) > 0. We say that (a, A) 
is bounded if sup„ >0 sup xeR / z \\A n (x)\\ < oo. 

Given two cocycles (a,A^) and (a,^4^ 2 ^), a (complex) conjugacy between them is a continuous 
B : R/Z — > SL(2, C) such that (| 1 . 3[) holds. The Lyapunov exponent is clearly invariant under 
conjugacies. 

We assume now that (a, A) is a real cocycle, that is, A £ C°(R/Z, SL(2, R)). The notion 
of real conjugacy (between real cocycles) is the same as before, except that we ask for B 6 
C°(R/Z, PSL(2, R)). Real conjugacies still preserve the Lyapunov exponent. An example where 
it is useful to allow B : R/Z -> PSL(2,R) instead of requiring B : R/Z -> SL(2,R) is given by the 
following well known result (see [MP] for the case of continuous time) . 

Theorem 2.1. Let (a, A) be a uniformly hyperbolic cocycle, with a Diophantine and A analytic. 
Then there exists an analytic B : R/Z — > PSL(2,R) such that B{x + a)A(x)B(x)^ 1 is constant. 

One can not always take B : R/Z -> SL(2,R) in Theorem [2~l1^ 



''Indeed, the continuous splitting of R 2 associated to a real uniformly hyperbolic cocycle may be topologically 
non-trivial (see §4.3 of [H] for an example). 
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We say that (a, A) is (analytically) reducible if it is (real) conjugate to a constant cocycle, and 
the conjugacy is analytic. 

The fundamental group of SL(2,R) is isomorphic to Z. Let 



(2.5) Re 



( cos 2tt9 — sin 2ir9 
lsin27r# cos27T0 



Any A : R/Z — > SL(2,R) is homotopic to x i— > R nx for some n G Z called the degree of A and 
denoted degA = n. 

Assume now that A : R/Z — > SL(2,R) is homotopic to the identity. Then there exists ip : 
R/Z x R/Z — > R and u : R/Z x R/Z ->■ R+ such that 

(2 - b > A{X) ysin 2th/ ) ~ U[x > v) ysin 2n{y + $(x, y)) 

The function ip is called a lift of A. Let /i be any probability on R/Z x R/Z which is invariant 
by the continuous map T : (x,y) i— » (x + a, y + ^>(a:, y)), projecting over Lebesgue measure on the 
first coordinate (for instance, take p as any accumulation point of i X)fe=o where v is Lebesgue 
measure on R/Z x R/Z). Then the number 



(2.7) p(a,A) = J ipdpmodZ 

does not depend on the choices of ip and p, and is called the fibered rotation number of (a, A), see 
[JM] and 0. 

It is immediate from the definition that 

(2.8) \p(a, B) - 9\ < C\\B - Rg\\ Q . 

The fibered rotation number is invariant under real conjugacies which are homotopic to the 
identity. In general, if {a,A^) and (a,A^) are real conjugate, B(x + a)A^ (x)B(x)^ 1 = A^\x), 
and B : R/Z — > PSL(2,R) has degree fc (that is, it is homotopic to i ^ Rkx/2) then 

(2.9) p(Q 1 4W) = p(Q,A( 2 ') + ifea/2. 

For uniformly hyperbolic cocycles there is the following well known-result (for a Diophantinc and 
A analytic, it is a consequence of Theorem 12. II) . see §5.17 of [H] , 

Theorem 2.2. Let (a, A) be a uniformly hyperbolic cocycle, with a G R\Q. Then 2p(a, A) G aZ+Z. 

2.2. Schrodinger operators. We consider now Schrodinger operators {H VyCC ,g}g £ s. ( we incorporate 
the coupling constant into v). The spectrum S = S v>ct does not depend on 9, and it is the set of E 
such that (ol,S v ^e) is n ot uniformly hyperbolic, with S Vj e given by (|1.4[) . 
The Lyapunov exponent is defined by L(E) = L(a, S Vi e)- 

Fixing a phase 9 and / G ^ 2 (Z), we let p = p{, a e be the spectral measure of H — B v . a fi 
corresponding to /. It is defined so that 



(2.10) ((H - E)- l f, f) = J ——dp(E') 

holds for E in the resolvent set C \ S. 

The integrated density of states is the function TV : R — ► [0, 1] defined by 

(2.11) N(E)= [ p f va9 (-^,E}d9, 

JM./Z 
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where / G ^ 2 (Z) is such that ||/||z2(z) = 1 (the definition is independent of the choice of /). It 
is a continuous non-decreasing surjective (for bounded potentials) function. The Thouless formula 
relates the Lyapunov exponent to the integrated density of states 

(2.12) L(E)= f In \E' - E\dN(E'). 

Jr 

There is also a relation to the fibered rotation number 

(2.13) N(E) = l-2p(a,S v , E ) 
where p(a, S ViE ) £ [0, 1/2]. 

2.3. The dual model. It turns out that the spectrum £ of H = H v , a ,e coincides with the spectrum 
£ of H = H v a g, a manifestation of Aubry duality (e.g., |GJLS| ). The spectral measures p, = fi v a 6 

can be defined using the analogous formula to (|2.10p . The integrated density of states N(E) = 
J p,g(—oo, E]dd then coincides with N(E). Berezanskii's theorem [Bel IS 1J gives in this context: 

Theorem 2.3. For any v, a, 9, f, and y? v a ^-almost every E, there exists a non-zero solution Hu = 
Eu with \uk\ < 1 + |fc|. 

Since S = £ is the union of the supports of the spectral measures, it follows: 

Theorem 2.4. For any v, a, 9, there exists a dense set of E G £ = T, ViCX such that there exists a 
non-zero solution Hu — Eu with \uk\ < 1 + |fc|. 

2.4. Localization and reducibility: dynamical formulation of Aubry duality. We describe 
the connection between localization and reducibility we mentioned in the introduction. The Theorem 
below is essentially a careful dynamical formulation of the classical Aubry duality. It appears in a 
similar form in P2 . 



Theorem 2.5. Let and let v : R/Z -> R be analytic. Let 9 G R, E G R be such that there 

exists a non-zero exponentially decaying solution Hit = Eu, and let A = S Vj e- 

(1) If 29<£a1 + Z, there exists B : R/Z -> SL(2, R) analytic such that B(x + a)A{x)B{x)- 1 = 
R±g (so (a, A) is reducible). In particular \\A n (x)\\ — 0(1), x G R/Z. 

(2) If29-kaeZ for some k G Z, there exists B : R/Z -> PSL(2,R) and k : R/Z -> R analytic 

such that B(x + a )A(x)B(x)- 1 = In Particular \\A n (x)\\ = 0( n), x G R/Z. 

If moreover a G DC, then n can be chosen to be a constant (and (a, A) is reducible). 
In either case p{a, S v> e) = ±$ + ma/2, for some m G Z. 

Proof. Let u{x) = Y,u k e 2mkx and U(x) = f e2 ™u( x )\ Then ^(a;) ■ U(x) = e 2 " e U(x + a). Let 
B(x) be the matrix with columns U(x) and U(x). By minimality ofim + a, det B is a constant. 

/ e 2iri0 q \ 

(Case A) If deti? ^ 0, we have B(x + a)^ 1 A(x)B{x) = -2-nie )■ ^ i s eas Y to see that 



e 

r ±^ 



det B = ±ci for some c > 0, we then take B = ,„ h ,„ B i , 

(2c) 1 / 2 M zp z 

(Case B) If det£? = then U(x) = ip(x)W(x) with W(x) a real vector defined up to sign and 
^(rr)! = 1. By minimality of x i— > x + a, W{x) ^ for every x G R/Z. The matrix B(x)^ 1 G 
PSL(2, R) with columns W and ^ xW R 1/4 W(x) is thus well defined. Then B(x + a)A(x)B(x)- 1 = 



(V d(x)-i)> where d{x) = ^m 1 ^ 6 - Since 



1 and d{x) is real, d(x) = ±1. If a G DC, 
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we can also further conjugate A to a constant parabolic (or identity) matrix by solving (using, say, 
Fourier series) the cohomological equation ±(f>(x + a) =p</>(x) = k(x) — J n(x)dx with J Q cj)(x)dx = 

in M/Z. Letting B'{x) = (\ B(x), we get B'(x + a)A(x)B' (x)- 1 = f* 1 So < x ) dx 

Assume that 29 £ aZ + Z. Then we can not be in case (B): indeed e 27Tl9 = ± ^.+a) i m P ues (using 
Fourier series) that ip(x) = e~ mkx (notice that if) is well defined only in R/2Z) and e 27Tl6 = ±e mka , 
that is 29 = fcamodZ. Thus we are in case (A), and the first statement follows immediately. 

Assume now that 29 — ka G Z. If we are in case (B), the second statement follows immediately. 
Otherwise we are in case (A), and we have B(x+a)A(x)B(x)~ 1 — R±g for some B : M./Z — > SL(2, R). 
We then set B'(x) — R-^k 2L B{x), and we get B' (x + a)A(x)B' (x)^ 1 = ±id. So the second statement 
still follows in this case (with k(x) = 0). 

The statement about the value of p follows immediately from (|2.9[) . □ 

If 29 € aZ + Z, we will say that 9 is rational (with respect to a). 

Remark 2.2. It is clear from the above proof that if the Fourier transform of u is analytic in a strip 
| Sfcs | < e then the matrix B given in item (1) is analytic in the same strip. As for item (2), it is still 
possible to define a conjugating matrix with a definite complex extension, but one must be more 
careful in its definition. 

2.5. Rational approximations. Let q n be the denominators of the approximants of a. We recall 
the basic properties: 

(2-14) Iknalk/z = i nf l|fca|k/z, 

l<k<q n + i—l 

(2.15) 1 > 9„+i||9„a|| R/z > 1/2. 

One aspect of the "good distribution" of orbits {x + ja} q ^ 1 which will find repeated use in this 
work is the following estimate. 

Lemma 2.6 (Lemma 9.7 of |AJ1| ). Let a e K \ Q, x e K and let < l < q n - 1 be such that 
| sin7r(a; + loa)\ is minimal. Then for some absolute constant C > 0. 

(2.16) -C\nq n < V] In | sin7r(a; + la)\ + {q n - 1) ln2 < Chxq n . 

1=0 

l^l 



3. General setup and statements of the main (localization and dynamical) 

estimates 

3.1. Almost localization for every 9. Let a G R, 9 G R, eo > 0. We say that k is an eo-resonance 
if \\29 - ka\\ R/z < e^l fe l £0 and \\29 - ka\\ m = min b1 <| fe | \\29 - ja\\ R/z . 

Remark 3.1. In particular, there always exists at least one resonance, 0. If a G DC(k,t), \\29 — 
fea|U/z < e~ |fe|£ ° implies \\26 - ka\\ n/z = min|j|<| fc | \\20 - ja|| R / z for k > C(k,t). 

We order the eo-resonances — no < \ni\ < |?t,2 < .... We say that 9 is eo-resonant if the set of 
resonances is infinite. If 9 is non-resonant, with the set of resonances {no, . . . , rij) we formally set 
Tij+i = oo. The Diophantine condition immediately implies exponential repulsion of resonances: 

Lemma 3.1. If a G DC(k,t), then |n 3 -+i| > c\\2d — tij a \\^% — ce ce °'" J ', where c = c(k,t,€q) > 0. 
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Remark 3.2. In case \\29 — ^W||r/z = 0, (so that 9 is rational with respect to a) we have an 
especially strong resonance at £o- In particular, 9 is non-resonant, as there will be no resonances n 
with \n\ > \£q\. 

A simple Borel-Cantelli argument shows that the set of non-resonant 9 has full Lebesgue measure. 

Definition 3.1. We say that the family {H v ,a,e}eeM exhibits almost localization if there exist 
Co > 0, C\ > 0, eo > 0, ei > 0, such that for every solution u of H v ^ a ,eu — Eu satisfying u n = 1 
and \uk\ < 1 + |fc|, and for every Co(l + \rij\) < k < C ~ 1 \rij + i\ we have \uk\ < C\e~ tlk (where the 
rij are the e -resonances of 9). 

Remark 3.3. It is clear from Theorem 12.31 that almost localization implies localization for non- 
resonant 9 (slowly growing generalized eigenfunctions can be always shifted and normalized to match 
the definition). 

Theorem 3.2. If v : R/Z — > K is analytic and \X\ < Ao(wj3 then {Hx v ,a,e}deR is almost localized 
for every a E DC. For v(x) = 2cos27ra:, we have Xq = 1. 

This theorem will be proved in S}5] 

3.2. Bounded eigenfunctions for every energy. The next result allows to pass from "every 9" 
statements to "every E" statements. 

Theorem 3.3. If E E E then there exists 9 E M. and a bounded solution of H va gu = Eu with 
Uq = 1 and \u n \ < 1. 

Proof. It is enough to show this for a dense set of E in the spectrum. By Theorem 12.41 we may 
assume that there is a generalized eigenfunction u! with subexponential growth for some phase. Fix 
N > 0, e > 0. Let k t be a sequence such that \u' k . = max|j|<(j + i)jv \u'A. It follows that there exists 
some (and indeed infinitely many) i with \u' k . +i | < (1 + |. By shifting the phase and rescaling, 
we obtain, for every e > 0, N > 0, some phase 9 N,e such that there exists an eigenfunction with 
Uq' £ — 1 and |u^ ,c | < 1 + e, \n\ < N. Passing to the limit as N — > oo and e — > 

we get the desired eigenfunction^ □ 

Remark 3.4. The result can be generalized to some classes of continuous ergodic operators. For 
dynamical systems, corresponding results appear in the work of Mahe. Namely, if a cocycle is not 
uniformly hyperbolic then there exists a vector that is never expanded, either in the future or in the 
past. 

3.3. Main dynamical estimates: quantitative Aubry duality. We fix a, v, X as in Theorem 
13.21 For every E E E^.a, let 9 — 9(E) be given by Theorem l3.3FI Let {nj} be the set of resonances 
of 9(E). Let A — S\ Vt E- In what follows, C is a large constant and c is a small constant, which are 
allowed to depend on v, A, a, but not on E or 9. 



' The proof actually gives a quantitative bound on the dependence of Ao on the analytic extension of v. More 
precisely, what is needed is that supic^!^ |Aof(^)| < erjt for some < e < 1, where Co and fco are absolute constants, 
see J5TT51 . J5TT91 . ( T5T2tTl . 

8 An alternative proof is the following. Let p/q be close to a, and let 0' G R. The spectrum of H' = H ViP / qt gj is 
close to the spectrum of H in the Hausdorff topology. Let E' be close to E and in the spectrum of H' . Then there 
is a non-zero periodic solution to H'u' = E'u' . Changing 8' to 8' + kp/q, we may assume that u' = 1 and \u' n \ < 1. 
Taking the limit p/q — » a, E' — > E, a limit 8 of 8' , and a pointwise limit u of u' , we get the statement. 

Notice that 8(E) is not necessarily uniquely defined. There could be uncountably many resonant 8, but at most 
finitely many non-resonant 8, corresponding to the same E. This does not concern our arguments, we just fix some 8. 
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Theorem 3.4. Fix some n = + l < oo and let N = |n.j+i|. Then there exists $ : R/Z — > SL(2, C) 
analytic with ||$|| cn -c < Cn c such that 

(a,, * (l+<1 ,™-. = (7 e -L.) + (*g $»), 

(3.2) ||<Zi|U-o, ll?3|| c „-o, HftlU-o < Ge- cAr 
and 

(3.3) UglU-o <c e - c »(Mi+»))- c . 

Theorem 3.5. fix some n = \nj\ + 1 < oo and Ze£ N = \nj + \\. Let L^ 1 = \\29 — njCtW^/z, and 
assume that < L^ 1 < c0 Then there exists W : R/Z — * SL(2, R) analytic such that | deg W\ < Cn, 
\\W\\ C < CL C and \\W(x + a)A{x)W{x)~ 1 - R T e\\c < Ce' cN . 

Remark 3.5. For N — oo, Theorem 13 . 5 1 gives a quantitative version of the first statement of Theorem 
1231 

We will prove Theorem [33] in $7] and Theorem 13. 51 in Sj8] All spectral results are consequences of 
those theorems combined with Theorem 13.21 

3.4. Outline of the rest of the paper. Almost reducibility (Theorem II. 4p . and the direct proof 
of Theorem 11.11 are immediate consequences of Theorem 13.51 Estimates related to modulus of 
continuity, including Theorem 1 1.61 are obtained from Theorem 13.41 All are more or less immediate. 
We will discuss those consequences in Sj4j 

The technical core of the paper is formed by the proofs of Theorems 13.21 13.41 and 13.51 The 
dynamical estimates build on preliminary estimates and ideas developed in $6l but are otherwise 
independent. Localization (jjSJ) and quantitative duality f? ^6l71 and [5]) are independent. Section [5] 
uses the machinery developed in [BJ1] and some shortcuts from [AJ1] , The techniques and ideas 
developed in § £|6l7i and [8] are new and do not use any ideas/methods from the existing literature. 

4. Easy spectral consequences of the main dynamical estimates 

4.1. Almost reducibility. We will show the following precise version of Theorem 1 1.41 

Theorem 4.1. Assume a 6 DC,0 < |A| < Xq(v). There exists c > (depending on X,v,a) with the 
following property. Let A — S\ V ,E- 

(1) If p(a, A) is c-resonant, there exists a sequence flW : R/Z -> SL(2,R) such that B^(x + 
a)A{x)B^ n \x)- 1 converges to a constant rotation uniformly in {\^sx\ < c}. 

(2) If p(a,A) is not c-resonant and 2p(a,A) (/ ctZ + Z then there exists B : R/Z — > SL(2,R), 
analytically extending to {|Sz| < c}, such that B{x + a)A(x)B(x) is a constant rotation. 

(3) If 2p(a,A) 6 aZ + Z then there exists B : R/Z -> PSL(2,R) analytic such that B(x + 
a)A(x)B(x) is a constant. 

Proof. If E £ S, by Theorem HH 2p(a, A) e aZ + Z, and by Theorem O (a, A) is reducible. 

Let E 6 S. If 6* is not eo- reson a n t, by Theorem 13.21 there is an exponentially decaying eigenfunc- 
tion. Theorem 12.51 thus applies, and the result holds in all cases (using Rcmark l2.2[) . 

Assume that 9 is eo-resonant (and thus 26 £ aZ + Z). Applying Theorem 13. 5[ we get for every j a 
matrix : R/Z -> SL(2,R) such that ®(x) = B^ (x + a)A(x)B^ (x)- 1 satisfies \\<f>(x) - R^ e \\c < 
Ce~ cN and | deg B»' | < C(\rij |+1) (where N — |nj+i|). To conclude, let us show that p is c-resonant. 



It is likely that the result holds assuming only L _1 > 0. 
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By (EU), \p(a, A) ±6 + degS^a| < Ce~ cN . Using dHJ) and Lemma Owe get for large j 

(4.1) \\2p(a,A) ± ( nj ±2degB^)a\\ m > \\^29±n j a\\ 1l/z - C e - cN > cN- c - C e - cN > 0, 

(4.2) \\2p(a, A) ± (n 3 ±2degB^)a\\ m/z < \\ T 29 ± n ja \\ M/z + Ce~ cN 

< e - e ol™jl _i_ Ce~ cN < e _c l" J 'le~ c ^" J±2dcsB(3) '- ) . 

For j large, f|4. 2|) and the Diophantinc condition imply that p(a,A) has a c-resonance at ^(nj ± 
2degfiW). If the set of c-resonances for p(a, ^4) would not be infinite, then we would have 
\\2p(a, A) ± (rij ± 2 degSW)a|| R/z = for large j. But this contradicts l|Oj) . □ 

Remark 4.1. This result is new even in the perturbative regime. Though Eliasson's scheme provides 
a sequence of approximate conjugacies to normal forms, for badly behaved p there is only control 
in a shrinking sequence of bands, as typical for KAM analysis. In fact it does not imply our strong 
definition of almost reducibility of cocycles except when it actually implies reducibility. 

Remark 4.2. The analytic conjugacy given in the third statement for J?6S can be shown to extend 
also to a definite strip. For E ^ E, we do not get any estimates on the analytic extension. 

We also state separately the following statement, already obtained as a part of the proof of 
Theorem 14.11 By Theorem 12. 5[ if 9(E) (as specified in Theorem 13.31 see also the footnote in 
Section l3~3f is not eo-resonant, then p(E) = ±9 + ka/2 for some k. There are reasons to believe that 
this does not hold for all E or flFl 

Instead, for resonant 9 we have the following statement. 

Theorem 4.2. Fix £eS and some 9(E). Under the conditions of Theorem \4-l\ and with the same 
c we have that if 9 is eo-resonant, then p is c-resonant. 

4.2. Open gaps for the almost Mathieu operator (direct proof). By Theorem 12.21 in the 
closure of a component of M \ E we must have N(E) — 1 — 2p G ctL + Z. The dry Ten Martini 
Problem is the conjecture that the converse holds in the case of the almost Mathieu operator. Since 
the integrated density of states is the same for H and H it doesn't matter whether to prove such 
statement for A or for A -1 . Theorem II .11 is thus a consequence of the following. 

Theorem 4.3. Let v(x) = 2cos27rir, and let < |A| < 1, a Diophantine. If E G E is such that 
N(E) G aZ + Z then E belongs to the boundary of a component o/R \ E. 

Proof. By [Plj . it is enough to show that (a, S\ Vy E) is reducible. By Theorem 14. 1\ if this was not 
the case then p(a, S\ Vt E) would be c-resonant. By the Diophantine condition and Remark 13. 2\ 
2p ^ aZ + Z. Since N = 1 — 2p, this is a contradiction. 

□ 

4.3. Precise bounds on growth, complex perturbations. Conventions below are as in section 
1331 

Theorem 4.4. Let n = \nA + 1< oo and let N = \n ]+1 \. There exists W = W (e) : M/Z -> SL(2,C) 
analytic such that, letting Z(x) — Z^(x) — W^(x + a)A(x)W^(x), we have 

(4.3) ||W|U-o < Ce- 1/4 and \\Z\\ cn -c < 1 + Ce~ cn(ln(n+1)rC e 1/2 , for Ce~ cN < e < crC c '. 

11 Otherwise the set of _B's corresponding to a given 9 would be countable even when Hg has singular continuous 
spectrum. If we relaxed slightly the definition of 9(E) to require |«o| > 1 — e instead of |«o| = 1 we could have argued 
that this is a contradiction since it can be shown in certain cases that all generalized eigenfunctions are bounded, and 
therefore for all energies E in the support of singular-continuous measure for a given 9 one can find 9(E) of the form 
9 + ka for some k. Such relaxed definition of 9(E) could also be used in our proofs, with only small changes. 
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Proof. Let <fr be given by TheoremEU Let D = ( ^ rf M where d = H^H^-ce 1 / 4 . Let W{x) 
D$(x). If e < cn' c we have ||W|| c „-c < Ce" 1 / 4 . Since 

(4.4) D ( ai a2 )D-i=(«} 
we get 

(4.5) Z( X ) = W(x + a)A( X )W(x)-* = g _L) + 7 ' 

with ||zi|| c „-c,||z 3 || c „-c,||z4|| c „-c < Ce-^e- cN and ||z 2 || c „-c < Ce^e-^M^+i))- . If £ > 
Ce~ cAr then ||Z|| c „-c < 1 + Ce-^^+^^e 1 / 2 . □ 



The following gives a direct proof of a perturbative estimate of [Ej (see Theorem IA.1I in the 
Appendix) . 

Corollary 4.5. For every s > we have \\A s \\q < C(l + s). Moreover, if 9 is non-rational then 

||i4.|| = o(l + 8). 

Proof. By Lemma FSTTT for every s > C, e — 1/s 2 is in the range specified in (|4.3[) for some choice of 
n = \rij\ + 1. Then 

(4.6) ln\\A s \\ < 21n||W||o + aln||Z|| < C - ~ hie + s c e - cn < ln (" +1 "' C e 1 / 2 <C + ln(l + s). 

This gives ||A S || < C(l + s). 

If 6> is non-rational and non-resonant then by Theorems 13.21 1231 ||A S || =0(1). If 9 is resonant 
then we have n — > oo as s — * oo, and we get the estimate 

(4.7) sup In \\AjW <C + ln(l + s), 

0<J<ce t! * , < ln ")~ C 's 

which implies the second statement. □ 

Corollary 4.6. If B : R/Z -> SL(2,C) is continuous then L(a, B) < C\\B - A\\l /2 . 

Proof. It is enough to consider the case when e = \\B — A\\o is sufficiently small. Then e is in 
the range specified by (|4.3[) for some n = \rij\ + 1. Let -B(x) = W(a; -I- a)i?(x)I4 7 (x)~ 1 . Then 
||S|| < ||Z|| + ||W||g||B - A|| < 1 + Ce 1 / 2 . Then L(a, B) = L{a, B) < In ||B|| < Ce 1 / 2 . □ 



As before, the estimate is improved when 9 is not rational. 



4.3.1. The integrated density of states: direct vroof of Theorem \l.b\ Theorem 1 1.61 follows easily from 
Corollary gH Indeed by the Thouless formula L(E) = j \n\E - E'\dN(E'). Thus L(E + ie) > 

L(E + ie) - L(E) = |/ln(l + m=WV ) dN ( E ') > c(N(E + e) - N{E - e)) for every e > 0. By 
Corollary HH L(E + ie) < Ce 1 / 2 , E E S. Thus N(E + e) - N(E - e) > Cc^e 1 / 2 for every 
0<e<l,£J€E. Since A" is locally constant in the complement of S, this means precisely that N 
is l/2-H61der. 
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5. Almost localization: proof of Theorem 13.21 

Although it will not be needed for the rest of the paper we will consider a weaker Diophantine 
condition on a. For v > 0, £ > 0, let EDC(V, £) be the set of a such that 

(5.1) \qa-p\>ve~ iq , pGZ,gGZ\{0}. 

Clearly for all k, t, £ > there exists v > such that DC(k, t) C EDC(V, £). If a £ EDC(^, £) then, 
by (|2.15j) . lnq n+ i < £q„ - Inu. 

We will prove the following precise version of Theorem 13.21 

Theorem 5.1. There exists Xo(v) > such that if < A < Ao, Co > 1, there exist £ = £(A, u, Co) > 
0, eo = eo(u,A) > 0, e\ — ei(w,A,Co) > 0, such that if a G EDC(^, £), t/iera i/ie family H\ VtQ , g 
is almost localized with parameters Co,eo,£i,Ci, where C\ = Ci(v, A, Co, v) > 0. For v{x) — 

±2 COS27T2, Ao = 1. 

5.1. Proof of Theorem 15.11 To simplify the notation, fix some v : R/Z — > M analytic, and set 
H = ff Q-e > = \H\ v a g. Through this section, C denotes an absolute large constant, while C„, 
for instance, denote a large constant that only depends on a. We warn that this is not the same 
convention used in the remaining of this paper. 

Fix an interval I CZ. Let F be the coupling operator between / and 1\I : 

Vi-j, Xi{i)+Xi{j) = 1) 
0, otherwise. 

Then we can write u = —{H — E — r) _1 ru from which for any x € I we obtain 

(5.2) u(x) = - ^ Gi(x,y)v y - k u{k), 

y£l,k£I 

where G/ is the Green's function of H restricted to the interval I : Gi — Gi(E) = (Ri(H — E)R* I )~ 1 
(where Rj : I 2 (Z) — > I 2 (I) and RJ : I 2 (I) — > / 2 (Z) denote the restriction and the inclusion operators). 
Set a k = ^ 

j:|i|>|fe|,jfc>0 

Fix E G M, m G M. A point y G Z will be called (m, fc)-regular if there exists an interval 
/ = [xx + 1, x% — 1], X2 = a;i + fc + 1, such that x\ < y < xi and 

^ |G J (i/,x)o x _ iei |<e-' nfc 

a;eJ,i=l,2 

Otherwise, y will be called (m, /c)-singular. 

The strategy is to show that the existence of a generalized eigenfunction as in Definition 13.11 
implies that k is (m,£(k)) regular for an appropriate m > and for £(k) comparable with k, with k 
in the desired "between the resonances" region. 

Define I C Z by I = [0, N — 1], N e N. We will omit the a, E dependence of various quantities in 
what follows, and all constants will be uniform for all E in the spectrum. We will also often assume 
v ^ 0, as otherwise our statements become trivial. Let Pn{&) = det Ri(H at g — E)R*j. Pn{0) is an 
even function of 9 + ^-^-ot and can be written as a polynomial of degree N in cos27r(# + ^-^-a) : 

Pn(0) = Ef =0 Cj cos^ 2tt(6 + £j±a) d = Qat(cos 2ir(0 + ^a)). 
Lemma 5.2. £ i In \P N (6)\d9 > - In A. 

Proof. It is proved by a standard Herman's subharmonicity argument [H] , D 
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Let A ktr = {9 e E, |Q fc (cos27r6')| < e (k+1)r }. The next lemma shows that every singular point 
"produces" a long piece of the trajectory of the rotation consisting of points belonging to an appro- 
priate Ak, r - It is rather immediate in the almost Mathieu case, and we will adapt the argument of 
|BJ1) for the general case. 

Lemma 5.3. There exists Xq — Xq(v) > such that for < A < Ao, < 6 < ^0 e = e(A, v) > 0, 
c = c(A, v) > 0, K C M \ Q compact, a 6 K, if y € Z is (c, N)-singular, N > N(X, v, K) and ieZ 
is such that y — (1 — S)N < x < y — JAf, we have that 9 + (x + N ~ 1 )a belongs to Ajy,— jnA— e- ^o 7 " 
v(a:) = ±2cos27ra, Ao = 1. 

Assume, without loss of generality, that Co\rij k \ < k < -^\nj k+ i\, and k is large (depending only 
on Co). We will define scales n > and s > 1 associated with k so that 

(5.3) 2sq n <(k< min{2(s + l)q n , 2q n+ i}, 

where £ = if 2|rij fc | < k < ^|n. Jfc +i| and £ = j^r- otherwise. Note that s,n depend on eo. 

Lemma 5.4. Assume there exists u as in Definition \S.1\ For A, c as in Lemma \5.S[ there exist 
£ = A, Co) > 0, 6q — 6q(v, A) > 0, smc/i that if a £ EDC(^, £), k > k(v,X,c,Co,v,£) and s, n 
are as above then k is (c, 6sg n — 1) -regular. 

Since sq n > Cq k the theorem now follows immediately from the definition of regularity and 
(15. 2[) . It therefore suffices to prove Lemmas 15-31 15.41 



5.1.1. Proof of Lemma \5.3\ . Without loss, set x = 0. Set x\ — —1, X2 — N. Assume 9 + ^j^-a 
Aiv,-inA-e! that is Pn(9) > X N e~ eN . We need to show that for y 6 [xijia] with dist(y, d[x\, X2]) > 
SN we have 



(5-4) (*)= \Gi(y,z)\\a z . Xi \<e 



cN 
zel, i=l,2 



By Cramer's rule Gi(y,z) = p^fo where \i y ^ z is the corresponding minor. The following lemma 
reduces the study of Green's function to the study of determinants of restrictions of H — E. 

Lemma 5.5 (Lemma 10, |BJlj ). 

M 

(5.5) [i y ^ = ^ a 7 det Ri\ 7 (H - E)R}^ |w 7i+1 - 7i | 

7 i— 1 

where the sum is taken over all ordered subsets 7 = (71, . . . ,7„) of L with 71 = y and 7„ = z, \"f \ — 
n — 1, and a 7 G { — 1,1}. 

The following is a simplified, yet more general, version of Lemma 11 of BJ1 . It gives an upper 
bound on [i VtZ . 

Lemma 5.6. For any Ac/, for sufficiently large N > N(v, X,K), 

(5.6) I deti? AA (H - E)R* AA \ < X- N e c ^ 2 ^N ^ + C^X^^j ^ . 

Remark 5.1. Unlike the corresponding upper bound in [BJlj this Lemma does not require a G DC. 



12 As will be clear from the proof, the result still holds for any < 8 < i, but the constants, including Ao, become 
dependent on S. 
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Proof. By Hadamard's bound, we have 

(5.7) \det R IXA (H-E)R* XA \< J[ ((2A" 1 cos 2tt(0 + 3 a) - E) 2 + \\v\\ 2 L 2 ) 1 ' 2 . 



jei\A 

Thus 



(5.8) m | det i?j\ A S).?^ | < ^ i ln((2A~ 1 cos 2n(6 + 3 a) - E) 2 + ||w||§) = ^ u(x + ja). 



2 

ie/\A jei\A 



Set At = {x : |2cos(27rx) — XE\ < t}, < t < 2 + [A-E7}, and let to be such that the Lebesque measure 
of A to is equal to Then by unique ergodicity of irrational rotations and continuity of u, we have 

(5.9) ^2 u(x+ja) <^2u{x+ja)- u(x+ja)+No(l) < N( u(x)(l-XA to )dx+o(l)), 

jei\A jei jei, 

x+ja£A t0 

Notice that a direct computation gives 



(5.10) / u{x)dx = In A" 1 + \ \n\z+ \J z 2 - 1||, 



where z = ^4t-t^ (see |BJlj for details). From this and \E\ < 2A 1 + ||i>||o in the spectrum, we 



get the estimate J u(x)dx < In A 1 + CA 1 / 2 ||w||q^ 2 . To complete the proof, it suffices to show that, 



letting s = &£, 



(5.11) / u(x)xA to (x)dx> S H\\v\\ + C- 1 X- 1 s 2 )-C\\v\\ 1 /2 X 1 ^. 



Let us split A t into four segments Ij of length Sj < 1 such that x 1— » (2cos(27ra) — XE) 2 is 
monotonic in Ij. Then, since 1 — cos27ra; > 16a; 2 for x < 1/4, 



(5.12) u{x)dx> -ln(|H| 2 , + 4A- 2 (l-cos(27ra)) 2 )dx> / -\n (\\v\\l + 256X- 2 x 4 ) dx 
Jij Jo 2 Jo 2 

> i S ,ln(|| V || 2 + 256A- 2 S |) -2s,. 
Since iniln (||w||o + 256A~ 2 ir 4 ) is concave on M + , we get, summing over j, 

(5.13) / u{x)dx>-sln(\\v\\ 2 + X- 2 s 4 )-2s. 

JA t0 2 

Considering separately the cases A _1 s 2 > c||i>||o and A _1 s 2 < c||t>|| we obtain that this implies 
(ISTTTD . □ 



Our assumption |P/v(6*)| > A - N e~ c , implies, with the notation of (|5.4[) and using (|5.5[) : 

N-l 



(5.14) (*)<(XeYJ2 E \ detR i\^H-E)Ri\^ Xi - lhl+1 \l[\ 



n=l i=l,2, 
7: |7|=7 



N - 1 / ^ + p2\-("+l) 

N 2 



71=1 2 = 1,2. ^ 



e -erb(7,i) 



7: |7|=n 

where < a < a(v) < 1 is such that 

(5.15) \v k \ <C7^e- |fc|CT , 



18 ARTUR AVILA AND SVETLANA JITOMIRSKAYA 

hi 



and 6(7, i) = |7| 7 |+i J2[=i l7i+i "Til- Let G M = {j, \l\ = n and 6(7,1) = b}. Then 

(5.16) (*) < e^H!^ (l^llo + c-^- 1 ^^- ) e-*#G 6 ,„ 



/ /" I 1 \2 \ — V"-"!"- 1 -/ /L 



n=l b 

N - 1 / , ,, , 'I Y2 x - ( " - J! 

e 

n=l x " ' b,G b , n jt® 

If G b! „ 7^ then < max{dist(y, dl), n + 1} < b < (n + 1)N < N 2 . By Stirling formula, setting 
b = rN, n + 1 = sb, we have ( b n ) < CrNe^ s > N where <p(s) = -sins - (1 - s) ln(l - s). Thus we 
can estimate 

(si—I \ — 1 \ —rsN 

- X r 2 s 2 ) e""VW rW 

0<S<1 

The desired exponential decay of (*) follows if 

(5.18) (**)= sup e + CIMiyV/ 2 + (ln2C' + lnC t , (T +lnA-21n<Ss-- + ^ ) <5s < 0. 

0<s<l V S S J 

This condition is satisfied (for appropriate e) if 

(5.19) c = iT^C^o-- 3 

is small. Indeed, using that ||i>||o < cc j-° ; and that 4>(s)/s < 1 — Ins, 

(5.20) (**)< e+ (edej 2 - cr + (in C + In c + 3 In ^ - |-) <5s. 

In case = 2cos27ra the result in the form we need, for any < A < 1 is stated in Lemma 9.2 
in [AJlj (uniformity on E was not claimed but is automatic from the argument). The simplification 
in this case is mainly due to the fact that, as can be seen by a direct computation (but also follows 
immediately from Lemma l5T5"j) . \/J> y ,z\ = \Py-i(0)PN-z(0 + za)\ and we have a uniform upper bound 
of the form In |Pfe(0)| < C e ,K + fc(hi A -1 + e). It is obtained by identifying Pk(0) with the upper-left 
coefficient of the fc-th iterate of the almost Mathieu cocycle (a, S\-i VjE ), and applying the formula 
L(a, Sx-i v e) = — hi A for < A < 1 and all E in the spectrum |BJ2] , □ 

5.1.2. Proof of Lemma \5.4\ We first recall some concepts introduced in |AJ1| . Sec. 9. 

Definition 5.1. We will say that the set {#1, . . . , 0k+i} is e-uniform if 



(5.21) max max 1 f 

zer-i,ili=i,-,fc+i ■ LX 



'z-cos27r^)| <gfce 



ze[-i,i]i=i,— ,k+i fj^ \cos2n6j — cos27r^)| 

e-uniformity (the smaller e the better) involves uniformity along with certain cumulative repulsion 
of ±0i(mod l)'s. 

Lemma 5.7. Let i > 0. If 6\, . . . , 6k+i £ ^4fc,-inA-e and k > fc(e, e) is sufficiently large, then 
{61, . . . , Ok+i} is not (e — 6)-uniform. 

Proof. Write polynomial Qk(z) in the Lagrange interpolation form using cos27r#i, . . . , cos27r0fe+i : 
(5.22) \Q k {z)\ 



S^nt 9 ^ n^(^-cos2^) 
g Qk (C ° S ^ ] n^-(cos2^-cos2^) 
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Let 9q be such that \Pk(6o)\ > A 1 . The lemma now follows immediately from (|5.22|) with z = 
cos27r(6»o + ^a). □ 

We now define I\ , 7 2 C Z as follows 

(1) Ji = [-2sq n + 1, 0] and J 2 = [fc - 2s 9n + 1, k + 2sq n ], if fc < l"*|±ii and > 0. 

(2) h = [l,2sq n ] and I 2 = [k - 2sq n + l,k + 2sq n ], if k < ^ 3 \ +l1 and n jk < 0. 

(3) h = [-2sq n + 1, 2sq n ] and J 2 = [* - 2sq n + l,k], if < fc < 



(4) A = [-2sg n + 1, 2sq n ] and J 2 = [k + 1, fe + 2sq n ], if fc > 



Kid 



2 



In either case, the total number of elements in 1\ U / 2 is 6sq n . Set 9j = 9 + ja, j 6 ii U J 2 . 
Lemma 5.8. If £ < C cq and k > k(Co, v,C) then set {9j}j^i 1 \ji 2 is Ceo + Cc £,-uniform. 
Proof. We will first estimate the numerator in (|5.21[) . We have, in each case, 

(5.23) In | cos 2na — cos 2-k9 2 -| 



jei 1 ui 2 



= ^2 ln|sin27r a - 3 | + ^ In | sin 2^ ° - 3 | + (6sq n - 1) In 2 

j'e-Tiu/2 je/iUJ 2 
= £+ + £_ + (6sg„ - l)ln2. 
Both E + and £_ consist of 6s terms of the form of (|2. 16[) plus 6s terms of the form 

(5.24) In min | sin(27r(x + — ))|, 

j=l,...,g n _i 2 

minus In | sin sM±\. Therefore, by (|2~T6|) 

(5.25) In | cos 2na — cos 2ir9j | < — 6sq n In 2 + Cs In q n . 

je-fiuj 2 

To estimate the denominator of (|5.2ip we represent it again in the form (|5.23[) with a = 9i. Then 

(5.26) £_ = ln l 

sin7r(i — j)a\. 

je/iu/ 2 

£_ consists of 6s terms of the form of (|2 . 16[) plus 6s — 1 minimum terms (since, when splited into 6s 
sums over intervals of length q n each, one of the sums will be exactly of the form of one of (|2.16p ) . For 
each i, j e I\ U I2, we have \i — j\ < k + Asq n < Cc sq n . In particular In | sin ir(i — j)a\ > —Cc sq n £,, 
for large q n (depending on v, £) . Notice that 

(5.27) max{ln | sinx|, In | sinir + ^ja\] > — 2^q n , x E K, < \j \ < q n +i, 
provided q n is sufficiently large (depending on v and £). Using that sq n < q n +i we get 

(5.28) £_ > -6sq n ln2-C Co sq n Z. 

Similarly, S + consists of 6s terms of the form of (|2.16[) plus 6s minimum terms, each of the form 

(5.29) In I sin(27r(0 + ^ + ^ a ))| 

for some |j| < Cc sq n < Cc q n +i, minus In | sin 2it(9 + ia)\ (which cancels a possible minimal term 
with i = j). Using (|5.27[) . we see that at most 6 minimum terms are smaller than —2^q n . Let us 
estimate the smallest term with j ^ i. 
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Consider first cases 1 and 2 of the definition of I\, 1^. Then, by the definition of h, i = f,2, 
and using (|5.3j) we have i + j =/= nj k and \i + j\ < \n,j k +i\. Therefore the smallest term is bounded 
below by —Ceosq n when k > 2\rij k \, and by —Cc„£,sq n , when k < 2\nj k \. In cases 3 and 4 we have 
i + j =/= rij k+ i, \i + j\ < 2\n,j k+ i\, and sq n > C^\nj k+ i\. Therefore the smallest term is bounded 
below by -C^\n jk+ i\ > -C Co £,sq n . 

Putting all together, we see that 

(5.30) £+ > -6sq n In 2 - Cs In q n - (Ce + C Co £)sg„. 

Combining, fOS]) . ([ST2^|) , and (jOOj) , we obtain 

max TT If Z COs27r ^l ^ e (gco+Cc b Q Jg »+C«lng n 

je/lU/ 2 I COS 27T0.- — COS 27T^ I 

as desired. □ 



We can now finish the proof of Lemma 15^41 By Lemmas 15 . 71 and 15.81 at least one of 9j, j G I\ U/2, 
is not in ^46sij„-i.-inA-(Cc 4+Ceo)-e where e can be made arbitrarily small for large n. From (15. 2p 
we see that the existence of a generalized eigenfunction u (as in Definition (|3.1|) ) implies that 
y = is (m, fc)-singular for any m > and for fc sufficiently large depending on m. Thus by 
Lemma [5.31 with sufficiently small eo, £ (depending also on Co), y = and <5 < |, we have that 
for all j G I\, Oj G A 6sqn _ 1 _ in x-{c Co ^+Ceo)-e f° r sufficiently large n. Let jo G ^2 be such that 
9 jo £ A 6sqn _x _i nX _(c Go t+Ce )-e- Since \j - k\ < 2sq n we can apply again Lemma [O] with S < ± 
to obtain that A: is (c, 6sq„ — l)-regular. □ 

6. Preliminary dynamical estimates 

6.1. Lagrange interpolation argument. Given Fourier coefficients w = {wk)k£Z and an interval 
I C Z, we let w 1 = Ylkel Wke 27Tikx . The length of the interval / = [a, 6] is \I\ = b - a. 

We will say that a trigonometrical polynomial p : R/Z — > C has essential degree at most k if its 
Fourier coefficients outside an interval / of length k are vanishing. 

Theorem 6.1. Let 1 < r < [q n+ \/ q n \. If p has essential degree at most k = rq n — 1 and xq G R/Z 
then for some absolute constant C 

(6-1) |b||o < Cq^ sup \p(x +ja)\. 

0<j<k 

Proof. We may assume that p(x) = P(e 27Ttx ) where P is a polynomial of degree k. Then by Lagrange 
interpolation, 



(6.2) p{x)=Y J P(x Q +ja) J] 



j=0 0<l<k. 



Thus 

(6.3) In |||)||o < lnrq„ + In sup \p(x+ja)\+ sup 2, m ' 



< , |1 _ p 2-!Ti(—ja+la) I ' 

o<j<fc o<j^fc, <l<k \i e 'I 

In order to prove (|6.ip it is enough to show that for l<s<r,0<j< rq n — 1 and x G R we 
have 

(6.4) ln\l-e 2 ™ (x+la) \ <Clng„, 

(s-l)q n <l<sq n -l. 
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(6.5) M 1 - e 27Tl{ - ja+la) \ > -Clnq n+1 . 

(s-l)q n <l<sq n -l, 

Consider first the case (s — l)q n < j < sq n — 1. Then 

(6.6) ln|l-e 27r ^ +/Q )| < ln I 1 ~ e 2 ^ x+la) \, 

[s-l)q n <l<sq n -l. (s-l)q„<l<sq n -l. 

where |1 — e 2m i x + l oa) | i s minimal. Using that |1 — e 27lly \ — 2| sin7ry| and Lemma \2.&\ we get (|6.4[) . 
The same argument gives (|6.5|) even more directly, with — Clnq n in the right hand side, since the 
sum to be estimated is already of the form considered in Lemma IzTlil 
Consider now the case s — 1 ^ [j/q n ]- Writing 

(6.7) Yl ln|l -e^^+'^l = ^ ln |1 - e 2 ™(*+ ; «)| + l n |1 - e ^i{x+i oa )^ 

(s-l)q„<l<sq n -l. (s-l)q n <l<sq n -l 

where (s — l)q n < lo < sq n — 1 is such that |1 — e 2iri(x+i a) | j g minimal, we see that (|6.4p follows 
from Lemma \2. 61 To obtain (16.51) from Lemma [2"T6l we must also show that 



(6.8) inf hill - e 27ri( '^' )Q | > -Cmg„+i. 

(s-l)q n <l<sq n -l 

But this follows from pTTlj) and (PT5)l . This concludes the proof of IpTf) . □ 
The Diophantine condition a 6 DC(k,t) implies 

(6.9) g„ +1 < k" 1 ^- 1 - 

In this case, (I6.1|) implies, with c = — > and C — C(k, t), 

(6.10) HpIIo < Ce c(1+fc)1_cin(1+fe) sup \p(x+ja)\, 

0<j<k 



-1„T-1 



since 1 + k = rq n < q n +i < K q\ 

6.2. Polynomial growth. From now on 

a 6 DC(k, t) and the family {H\ v ,a,e}ees. is almost localized with parameters (eo, Co, Ci, ei). 
Let -E be in the spectrum and let A = S\ v ,e- Let 9 — 0(E) be given by Theorem 13. 3[ and let {rij} 
be the set of resonances of 9(E). 

All constants may depend on k, t, Co, eo, ei, Ci, and on bounds on the analytic extension of Xv. 
We will use C to denote large constants and c to denote small constants. Further dependence on 
other parameters, will be explicitly indicated, for instance, we will use C$ for a large constant that 
depends on all parameters above and also on an arbitrary parameter 5 > 0. 

The specific values of generic constants such as C, c, Cg, may change through the argument, even 
when they appear in the same formula. In a few places we will use non-generic constants, denoted 
C", C^> to simplify later referencing. 

Theorem 6.2. We have L(a, A) = 0. 

Proof. Fixing a non-resonant 9' , Theorem 13.21 fsee Remark |3.3[) implies localization for H\ V: u,e>- By 
Thcorem l2.51 L(a, S\ Vj e>) = for a dense set of E' in the spectrum. By continuity of the Lyapunov 
exponent, |B J2] . this implies that L(a, S\ v- e) = for all E in the spectrum^ □ 



^An alternative argument for zero Lyapunov exponent to be dense is the following. By duality IGJLSI . if H has 
pure point spectrum for almost every then H has absolutely continuous spectrum for almost every 6 (see discussion 
in | BJ1I for the models considered here). By the Ishii-Pastur Theorem (e.g.. [CLl ). the Lyapunov exponent is zero 
densely in the spectrum. 
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In our context (see Remark l2.1j) . this means that 

(6.11) sup ||4t(a:)|| < C s e Sk , S > 0. 

£ceR/z 

Our goal in this section is to improve subexponential growth to polynomial growth. This is of 
course still not optimal (see Corollary 14. 5p . but it is a good starting point. 

Theorem 6.3. We have \\A S \\ < C(l + s) c , s > 0. 

We note that the proof of this estimate involves some themes which will appear again later. 
Choose 4Cb(|nj|+l) < m < |nj+i| of the form m = rq k ~l < qk+i, let / = [— [m/2], m— [m/2]] 

and define u(x) — u T (x). Let U(x) = I , \ )• Then 



2ni0TT(„ i „A 47ri9 I h(x) 



u(x — a) 

(6.12) A{x) ■ U(x) -e 27rl6 U{x + a) = e , (| 
where 

(6.13) h k = x/(fc)2cos27r(6» + ka)u k + ^xi(k - j)vjUk-j, 
where xi is the characteristic function of /. Since Hu = Eu, we also have 

(6.14) -hk = Xz\/(fc)2cos27r(6! + ka)u k + Xz\i{ k ~ j)vjUk-j- 

The estimates \v,k\ < Cie~ ei ' fc for m/4 < |fc| < m, \u k \ < 1 for all k and \v k \ < Ce~ c l fc l for all k 
then imply that \h k \ < Ce- cm e~ c ^, that is 



-cm 



(6.15) ||ft|| c <C'e- 
Theorem 6.4. We have inf xeK /z \\U{x)\\ > c s e- 2Sm , 5 > 0, m > C. 

Proof. Otherwise, by (jlTTT]) . (j!Tl"2")) . (|6"T5|) . \u(x + ia)\ < e"" 5 ™ for some x € M and < i < m. Then, 
by TheoremO ||uj| < Cq^ r e~ Sm < C s e~ 5m / 2 . This contradicts J u(x) = 1. □ 

Let B(x) be the matrix 

'e 2 ™ e u{x) - w l xW u(x - aY 



(6.16) 




By Theorem 16.41 and the trivial estimate \\U\\o < Cm we have 
(6.17) ||fl|| < C s e 2Srn , 5>0. 

Thus 



e 2 " 9 \ fpiix) b(x) 
e" M + UW Hx) 



(6.18) B(x + a)" 1 ^)^^) 

with 

< C 5 e 45 ' 



(6.19) 



/3i(x) b(x) 
fo{x) p 4 (x) 





Since the first column of B satisfies (|6.12|) and -B £ SL(2, C), ()6.15|) and (|6.19|1 for appropriate 6 give 
(6.20) ||/3i||o,||/?3||o,||/34||o<C e - cm . 
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Taking $(a;) the product of B(x) and a constant diagonal matrix, = DB(x) , where 

D = (n w-ij ' witn rf2 = max {ll/33||o /2 , e~ m } we get 



d~ 

2rrie 



/ 2-KiB n \ 

(6.21) *{x + aJA^)*^)" 1 = ^ Q e _L e J + Q(a), 

where, by (jOU]) . (|637jl .([6TTg ]) .and (jOll with appropriate 6, \\Q\\ < Ce- cm and ||$|| < Ce m . This 
implies 



(6.22) ||A fl || < 2||$||S < C7e° m , < s < c w e c m . 

Proof of Theorem \6.3\ Let c' 1 ' be as in (|6 . 22f) . For fixed s, let m be minimal such that s < <P^e c >m , 
4Co(l + | | ) < m < Cq and m = rq k — 1 < Qk+i for some j,k,r. By Lemma 13.11 

m< C* + C71n(l + s). By (ESI, ||A s ||o < C(l + s) c . □ 

6.3. Improved estimate on almost invariant sections. In the previous argument, polynomial 
growth was obtained from an estimate on a suitably chosen "approximate Bloch wave" . The first 
application of polynomial growth is that it allows us to deal, through a bootstrap argument, with 
much finer approximate Bloch waves. 

Lemma 6.5. Let n — \n,j\ + 1 < oo and let N = Let u(x) — u 1 (x) for I = \— Cq 1 N + 

_i fe 27Tie u(x)\ 
1,C N — 1J. Define as before U(x) — , {)■ Then we have the estimates 

27ri0rr/ , „.\ , ( H x )\ „„•,>, IIj.II ^ -cN 



(6.23) A[x) ■ U{x) = e 27Tl0 U(x + a) + \"y J , with \\h\\ c < e 

(6.24) inf \\U{x)\\>cn- c , N>C. 

\Qz\<cn c 

(6.25) ||t/|| < Cn 

Proof. Estimate (|6~23ll is obtained exactly as (|6T5]l . For (|6?24| . let 4C n < m < C^N be of the 
form m = q k - 1, and let J = [-[m/2],m - [m/2]]. Define f7 J (a;) = (1 J(/_ ^\ ) • Theorem HU 
implies 

(6.26) \\U-U J \\ C < Ce- cm . 
Arguing as in Theorem 16.41 we get 

(6.27) inf \\U J (x)\\ > cm- c , m>C 

(the estimate is better since |J| = q k — 1, and because we can use Theorem 16.31 instead of (|6.1ip ). 
Obvious bounds on the derivative of U J then give \\-^U J (z)\\ cm -i < Cm 2 , which together with 
(|6T2l)|) and fOT) imply 

(6.28) inf \\U(z)\\>cm- c , m > C (1) . 

\^sz\<cm~ c 

By Lemma |3. II and the Diophantine condition, if m > max{4Co^, C^} is minimal of the form 
m = q k -l then m < Cn c < C ~ 1 N, N > C. Together with f!T2"8"j) . this gives (j!T2"4"|) . 

Estimate (|6.25p is immediate by Theorem 13.21 □ 
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7. (Complex) almost triangularization: proof of Theorem 13.41 
Let U(x) be as in Lemma [6751 Let B(x) be denned by (|6. 16|) . By Lemma T6. 51 we get the bound 
(7.1) \\B\\ cn -o < Cn c . 

Here the complex extension of B is the holomorphic one (not the one given by (|6.16[) ). By the same 
computation as in (|6.18[) and using Lemma 16.51 we get the improved estimate 

(7-2) B(x + a) >A(x)B(x) = ^ Q ^ j + j*^ ^ j 

where ||6|| c „-c < Cn c and ||/3i|| C n~ c > 1 1 /?3 1 1 en— c » IIAtHcn- < Ce~ cN . Thus the Fourier coefficients 
of b satisfy the estimate 

(7.3) \b k \ < Cn c e- cn ~ Cw . 

Notice that the estimates so far cover the case of the first resonance (n = 1), so we will assume 
from now on that n > 1. 

In the estimates to follow, we shift one scale back. So take n = 1 + |n 3 -_i | , jV = |rij-|, and the 
remaining notation is clear. Split b(x) = b l (x) +b h (x) where the first term has the (low frequency) 
Fourier coefficients with \k\ < N_ and b h is the rest. We solve exactly 

/ 1 -^x + a)\( e 2 ^ 9 b\x)\fl <f>(x)\ _ fe 2 * ie 
{iA) \i) 1 ) \ e- 2 ™ e J\0 1 ) ~ { e- 2ml 

This corresponds to solving 

(7.5) b\x) - e- 2m9 ^[x + a) + e 2m \{x) = 0, 

or in terms of Fourier coefficients, 



for \k\ < N_ and <£ fe = for |fc| > JV. By (fT3jl . and since a £ DC(k,t) and |A;| < JV, we have that 
U{x)\\ C n-c < CN C . We conclude that, with $(x) = ^ Bix)- 1 we have 

(7.7) $( I + a)W)- 1 = ( e " e -L)+£W 



where 

and |||>]] c „-c < CJV C . Since, by (fT3|) . H&^U-c < C*n c e" c ^ C ^ < Ce"^^, we also have 
||£||cn-o < C*e- c ^ _C ^. It follows that 

(7.9) \\A s \\ CR -c <CN C , 0<s<ce^~ C ^. 
Back to the original scale, we get using Lemma |3~T1 

(7.10) \\A s \\ c(lnn) -c < Cn c , < s < ce cn ^ C . 

Let A > n be minimal so that |Sfc| < Ae~ A 1*1, k G Z. Let us consider a different decomposition 
6 = 6 r + 6 ; + b h so that now V has only the (resonant) n^-th Fourier coefficient, b l has the (low 
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frequency) Fourier coefficients with \k\ < A 3 (except for nj) and b h has the (high frequency) Fourier 
coefficients with \k\ > A 3 . By (|7.3[) . A < Cn c , and the definition of A gives 

(7.H) H^IL-c < Ce~ n . 

We can solve 



(7.12) 



1 -<j>(x + a)\ (e 2m9 b l {x) \ (l <f>(x)\ _ (e 2 ™ 9 
1 Ho e- 2 ™ 9 [O 1 j ~ \ e- 2 ™ 9 



As in {Li]) it follows that \\(f>(x)\\ cn -c < Cn c . We conclude that, with = K ^M -1 



__ , e 27Tie \ t fO b r {x) + b h {x) 



we obtain 

(7.13) <S>(x + a)A(x)<P(x)-' = P Q g _^ e j + I J " ^ ' Q " ^ 1 ) + *(,-) 
with 

(7.14) ||*IU-° < Cn c and ||*|| cn -o < C7 e - cAr . 
We now estimate b r {x) — b nj e 2 ^ lnjX . Let 



2-rrie 



(7-15) W(*) = (J Q e lH) 

We can compute exactly 



(7.16) W s 



Q ^ — I-KISQ 



with \b°(x)\ = \b nj j:t-le- 2 ^ 29 -^\ = \k /Zl S $t%$ \> if sinn(2e -nja) ? and \V(x)\ 
s\b nj \ otherwise. Therefore we have 



s\b nj \ 

'a||0 ^ 

On the other hand, by (|7.3p . 



(7.17) || W ,|| > 0<s<pe-n ja \\-} z /10. 



(7.18) ||W s || <l + s|6„ 3 .| <C*(l + s)n c , s > 0. 
Using ([7715]) . we get A = + a) _1 (VF(a;) + Z(z))$(x), so that 

(7.19) K||o > ||$|lo 2 Ml^llo-E (j)ll Z llo ^JI^IIo + ") ' 

By (|77l4"|) and (|7TT|) we have ||Z|| < Ce~ n . Thus, (fTTTEj) and (17713)) imply 

(7.20) ||A s ||o>cn- c '(||VF s ||o-C , e- cn ), < s < ce cn . 

By (177101) and ifT^Oj) . ||W a || < CVt c , < s < ce c "( ln ™r c . Together with Lemma O and (177171) . 
we get the estimate 

(7.21) \b nj \ < Cn c ' e - c "( ln »r c < c e - C «(inn)- c ^ 

The result follows from (|7TT3l) . ([77141) . ([TTTTj) and ([772T]) . □ 
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8. (Real) almost conjugacy to rotations: proof of Theorem 13.5 



Let U{x) be as in Lemma [6751 Let B(x) be the matrix with columns U(x) and U{x) on E/Z. Let 
= 1 1 26* — rija\\K/z- As standing hypothesis below, we assume 

(8.1) 0<L- 1 <c. 
Lemma 8.1. For any e > 0, we have 

(8.2) inf |detE(x)| > clT c , 

ifER/Z 

Proof. Recall that for any 2x2 complex matrix M with columns V and W, 

(8.3) |detM| = [[V||min||Vr-AV||. 

Notice that the minimizing A satisfies |[AV|| < ||W||. 

Minimize over A e C, x e E/Z the quantity \\e-* in * x U{x) - \e mn i x U(x)\\. This gives some A , 
Xq. If the estimate does not hold then, using Lemma [6.5[ Theorem 16.31 (|8 . 3[) . and Lemma f3.ll we 
would have 

(8.4) \\e- 2 ™ 9 e-™i X0 U(x~a + sa) ~ e 2ms6 'e m ' x °X U{x + sa)\\ < L~ l , < s < L. 



This implies that ||e-' ri7l ^ a! °+ aot )Z7( :c o + sa) - A e 7rm ^ :Eo + sa )[/(a; + sa)\\ < L^ 1 + CsL^WUWo < 
L- 1 / 3 for < s < L 1 / 3 (recalling ([BT25])). 

As in the proof of Lemma l6.51 let m > (InL) 2 > ACgn (here the second inequality is a consequence 
of (|8.ip ) be minimal of the form m = qk — 1, let J = [— [m/2],m — [to/2]] and let U J (x) = 

j, \ )• -By Lemma [3.11 and the Diophantine condition, m < C(\nL) c < C N. Theorem 

then gives \\U — U J \\ C < L^ 1 . Using Theorem 16.11 and Theorem 16.31 we obtain \\e~ 7rmjX U J (x) — 
e ™i*\ U J (x)\\o < CQjiL) c (L-^ 3 + 2L- 1 ). Thus 

(8.5) sup \\e-™* x U(xj -e^ x \ U{x)\\ < iT 1 / 4 . 

xGR/Z 



Substituting x x = x +sain JHU and taking s = [L/2], we get \\ie- 1Tin ^ Xl U{x 1 )+ie mn ' Xl \ U{x 1 )\\ < 
L- 1 + CL-^UWo, so that, by (JHISJ) and (fQ5l) . ||t/(xi)|| < i" 1 / 5 . 

However, (|6.24[) and Lemma EUJ imply ||L/(xi)|| > c(\nL)~ c , giving a contradiction. 

□ 



Lemma 8.2. Let x e E/Z. Tften 



(8.6) sup \detB(z) -det5(x ) < Ce 

Oz|<c 



cJV 



Proof. By A(x)B(x) = . " 9 B(x+a)+ ^ , ||A|| e < e - cAr , thus | det£(x - 



e -2^e)B(x+a)+i V Q ' 

a) - detfl(x )| < Ce- cN . This gives |detB(x + ka) - det£?(x )| < C e - cN ', < fc < 4iV. The 
function x i— > det £?(xo + fea) — det -B(x) is a trigonometrical polynomial of essential degree bounded 
by 47V. By TheoremO \detB(x) - detB(x )| < Ce~ cN ', x e E/Z. By Theorem[37J ||E/(z)|| c < 
Ce Cn so z ^ det B(z) - det£(x ) is bounded by Ce Cn over \$sz\ < c (2) . By the Hadamard 
three-circle theorem, 

(8.7) In sup \ det B(z) -detS(x )| < (1 -S)ln sup +Shi sup < -cN, < 6 < c. 

|Gz|=<5c< 2 > |9a|=Q |Gz|=c< 2 > 



14 A ctually, even a trivial bound like ||{/(z)|| c < Ce CN would do here. 
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Thus \det B(z) -detB(x )\ < Ce- cN , \$$z\ < c. □ 
Theorem 8.3. We have 

(8.8) inf |3fdet£(z)| > clr c > 

9z|<c A<° 

Proof. Notice that 'StdetB(x) = 0, x E R/Z. The result now follows from (jO]) and ([5TB]) (taking 
into account Lemma 13. II). □ 



Remark 8.1. Optimizing the method (and using other estimates obtained in this paper) gives 
| det-B(z)| > cgL^ 1 ^ 6 , |9fz| < c. This will be explored and used in |AJ2] , 

Take now S = WJ , T = on R/Z. Then B = (S ±T) ^ . _\j . Let W 1 be the matrix with 

/ 1 1 \ /e 27 " 9 \ / 1 1 \ _1 
columns S and ±T, so to have detWi > 0. Since I ^ J I ^ e -2ivi6 J ( _y ^ J = ^t#' 

by (|6.23[) we have 

(8.9) ||A(a;) • Wi(ic) - VKi(x + a) -i? Te || < C*e" cAr , < c 

(the complex extension considered here is the holomorphic one) . LetW(x) = | det-B(x)/2| _1 / 2 Wi(x) 
on R/Z, so to have det W = 1 (by Theorem l8.3l there is no problem with branching when extending 
|detB(x)| -1 / 2 to |3x| < c). Then 



(8.10) 



a, \ tt r i \ det£(x + a)| 1/2 T . , n 



< C*e~ cAr , |Qfa:| < c. 



A combination of (|8 - 6(1 and Theorem 18.31 gives 

detB(x + a)\^ 2 



.11 



< C7e" cAr , |3?x| < c. 



|det5(x)|V2 

Also, using Theorems O and [HJJ we get ||W|| C < CN c e cn . Thus 

(8.12) ||A(x) • W(x) - W(x + a) ■ R T g\\ < Ce~ cN , |3x| < c. 

To conclude, we need to show that |degVy| < Cn. Set U{x) = e m ^ x U{x), and S = W, f 
3*7, Wx = (S ±f), W(x) = \detB(x)/2\ 1 / 2 W 1 (x) on K/2Z. Then W(x) = W(x)R 

deg W-2 deg W = Tnj, where deg W is the degree of W : R/2Z -» SL(2, R) and deg W is understood 
in the usual sense, as the degree of W : R/Z — > SL(2, R). 

The degree of : R/2Z — > SL(2,R) is the same as the degree of any of its columns, considered 
as maps R/2Z -> R 2 \ {0}. Thus we need to show that the degree of M : R/2Z -» R 2 \ {0}, for 
either M = 5 or for M = T satisfies | degM| < Cn. 

Notice that 



t , ^ , so 

! 2 



(8.13) 11/ e- 7rm ^(S*(a;) +iT(x))dx\\ = || / [7(x)dx|| > 2, 

JR/2Z JR/2Z 

since u = 1. Select M(x) = S{x) or M(x) = f(x) so that / R/2Z ||M(x)|| > 1. Then, by §3^, for 

|9fx| < c we have \\A(x) ■ U(x) - e^ 29 - n ' a ^U(x + a)\\ < Ce' cN . By Lemma EQ and using (jOS]) 
and L- 1 < e~ eo l™l, we have \\\A(x) ■ M(x)\\ - ||M(x + a)||| < CL^n < CL~ C , x G R/2Z. Using 
Theorem 16.31 we obtain, arguing as in the proof of Lemma 16.51 



(8.14) inf \\M(x)\\>cn~ c . 

| <sx\ <cn~ ° 



2X 
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This already gives the bound |degM| < Cn by the obvious derivative estimate. To obtain a 
linear bound, let M(x) be a vector obtained by cutting off the Fourier modes of M(x) with ||| > A, 
where A is chosen minimal so that \\M{x) - M{x)\\ < ||M(a;)[[, x G R/2Z. Then by Theorem EH 
and (|8.14j) . A < Cn. By Rouche's Theorem, the degree of M is the same as the degree of M. 
Consider now a coordinate of x i— > M(2x) which is not identically vanishing. It is a trigonometric 
polynomial of essential degree at most Cn, and since it not identically vanishing, it has at most Cn 
zeroes in R/Z. It follows that | degM| (and thus | degM|) is bounded by Cn, and we conclude that 
\degW\<Cn. □ 

Appendix A. The perturbative theory of Eliasson 

The following result is due to Hakan Eliasson. 

Theorem A.l. Let a G DC(k, r), and let A : K/Z — > SL(2, R) be analytic. Assume that \\A— A*|| e < 
C ' 1 e ro where A* G SL(2,R), < e < 1, r Q = r (r) and C = C {k,t, ||A*||). TTien 

(1) // p(a,A) is either Diophantine or rational with respect to a then {a, A) is analytically 
reducible, 

(2) // (a, A) is not hyperbolic then sup r||Afc||o < Co, 

(3) If p{a,A) is not rational then lim ^||Afc||o = 0. 

Actually what is considered in [E] is the case of continuous time, and for cocycles of Schrodinger 
type (he also treats the case of several frequencies) . The considerations for discrete time are similar 
and are carried out in the thesis of his student Sana Ben Hadj Amor [A"m] . In [EJ, Eliasson goes on 
to establish absolutely continuous spectrum for all phases for the associated Schrodinger operators, 
while the proof of l/2-H61der continuity of the integrated density of states is the main result of 
Amor's thesis. 

The KAM scheme of Eliasson has been extended to the C°° case |AK2] . This generalization 
motivated the introduction of the concept of "almost reducibility" in the C°° case. 

A.l. Quick reduction to the perturbative regime. Reduction to the perturbative regime of 
Eliasson is a much less subtle result than showing almost reducibility. It can be concluded from 
almost localization just after fjHl since is clearlj0 implied by the following result. 

Theorem A. 2. If is resonant, there exists arbitrarily small e > 0, —1 < K < 1, and B : R/Z — > 
PSL(2,R) analytic such that \\B\\ t < Ce~ c and such that \\B{x + a)A(x)B{x)- 1 -A*|| e < Ce~ cj(C , 

where A jf ={\ ^ 

Proof. The beginning of the proof coincides with a small part of the argument used in the proof of 
Theorem 13 . 5 1 but we briefly repeat it here for the readers' convenience. Let n = \nj \ + 1, N = \rij+i\, 

and let I = l-C^N + l.C^N - ll. Let u = u 1 , and define U(x) = U ( X A as usua i. Let 

u u \u(x — a) J 

U(x) = e m i x ll(x) and let S(x) = m(x), f{x) = QU(x). Since (I5T31 holds, we can choose M = S 

or M = T so that / R / 2Z ||M(a;)|| > 1. Notice that by Lemma [6751 and the definition of resonance 

(A.l) A{x) ■ M{x) = M(x + a)+ 0(e~ c "), \%x\ < c. 

Then, arguing as in Lemma 16.51 we have (|8.14p . 



^^Recall that if 9 is non-resonant then by Theorems 12.51 and 13.21 (a. A) is conjugate to a constant cocycle (which 
is obviously in Eliasson's perturbative regime). 
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Let W be the matrix with columns M and j^^Ri/^M . Then by (|A.1|) and (|8. 14)1 we have 
W(x + a)- l A(x)W(x) = (l f) + &0*)j with \\p x \\ cn _ 0y \\p 3 \\ cn ^, \\04 m -o < C7e" c " 

and ||/?2 ||cri- c < Cn c . Solve the cohomological equation (j>(x + a) + /?2 (x) — <t>{x) = (32{x)dx/2 = b, 
with <p(x)dx = 0. Then, by writing out the Fourier coefficients and since a & DC we get ||</>|| crl - c 5= 

Cn c . Let $(x) = f J W(x)-\ then *(<b + aJA^a:)- 1 = (l + H where II- 3 < 

Ce- C ™. Let d 2 = min^l -1 , 1} and let D = f ^ . Then = + a)A(x)$(a;)- 1 £)- 1 is 

such that 

(A.2) ll-B-fii ^ IU-o < Ce" c ", 



for some — 1 < ft < 1. □ 
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